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The ABC model is a simple diffusive one-dimensional non-equilibrium system which exhibits a phase tran- 
sition. Here we show that the cumulants of the currents of particles through the system become singular near 
the phase transition. At the transition, they exhibit an anomalous dependence on the system size (an anomalous 
Fourier's law). An effective theory for the dynamics of the single mode which becomes unstable at the transition 
allows one to predict this anomalous scaling. 



PACS numbers: 02.50.-r, 05.40.-a, 05.70 Ln, 82.20-w 

A lot of work has been devoted recently to the study of the 
fluctuations of the current of heat or ofparticles through non- 
equilibrium one dimensional systems ill 4l2fl . In such studies 
the basic quantity one considers is the total flux Q{t) of en- 
ergy or of particles through a section of the system during 
time t. In the steady state this flux Q{t) fluctuates due to the 
randomness of the initial condition for purely deterministic 
models and due to the noisy dynamics in stochastic models 
(here we only discuss classical systems: see IllSl - Uql for the 
quantum case). If one assumes that the total energy or the 
total number of particles in the system remains bounded, the 
average current limj_j.oo \ as well as the higher cumu- 
lants limt_j.oo ^^^*^ '" of the flux Q{t) do not depend on the 
section of the system where this flux is measured. 

For a one dimensional system of length L, a central ques- 
tion is the size dependence of these cumulants lUTIl . In partic- 
ular one would like to know whether a given system satisfies 
Fourier's law, meaning that, for large L, the average current 
scales like 1/L; 



limi« 



L 



(1) 



where the prefactor Ai depends on the temperatures Ti and 
T2 of the two heat baths or on the chemical potentials jii and 
/i2 of the two reservoirs of particles at the ends of the system. 
At equilibrium (Ti = T2 or /ii = ^2) the prefactor Ai in ^ 
vanishes but the question of the validity of Fourier's law re- 
mains. One then wants to know whether the second cumulant 
of Q{t) scales like 1/L. 
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One can show that ^ holds for diffusive systems such as 
the SSEP (symmetric simple exclusion process)lill fl [isl [l9ll 
or the KMP (Kipnis-Marchioro-Presutti) modeling). The 
macroscopic fluctuation theory developed by Bertini et al. 
121, IZll I23I allows one also to determine |1| all the cumulants 
of the flux Q{t), with the result that they all scale with system 
size as \/L. 



Even corrections of order 1/L^ have been computed in some 
casesSH. 

For mechanical systems with deterministic dynamics, in 
particular systems which conserve momentum, the average 
current scales as a non-integer power of the system size: 

limM!))^4^ 
t^oo t L1-" 

The exponent a takes the value 1/2 for some exactly sol- 
uble special modelsyfl. Values ranging from 0.25 to 0.4 
have also been reported in simulations depending on the 
model considered II23I - I27II. T heoretical predictions based on a 
mode coupling approach ll28l 12911 or on renormalization group 
calculations JSOQ confirm this anomalous Fourier's law. Less is 
known on the size dependence of the higher cumulants, which 
are numerically harder to measure, except that they vary as 
power laws of the system size, with exponents which seem to 
depend on the geometry yH]. 

Here we consider the ABC model ll32l l33ll . a diffusive sys- 
tem which is known to exhibit a phase transition Il34l438ll : 
we study the fluctuations of the cuiTent near this transition. 
Generically, outside the transition the cumulants have a diffu- 
sive scaling ^. Here, we show that the amplitudes An be- 
come singular as one approaches the transition, and that the 
cumulants of Q{t) exhibit anomalous scalings at the transi- 
tion. When the transition is second order, due to the destabi- 
lization of a single Fourier mode of the density Il34ll , the fluctu- 
ations in the whole critical regime can be understood in terms 
of a Langevin equation for a single complex variable which 
represents the amplitude and the phase of this Fourier mode. 



I. DEFINITION 
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The ABC modeI ll32l438ll is a one-dimensional lattice gas, 
where each site is occupied by one of three types of particles, 
A, B and C. Neighboring sites exchange particles at the rates 

AB^BA 
1 

BC^CB 
1 

CA^AC 
1 



with an asymmetry g < 1 . Here, we consider the model on 
a ring of L sites. Since the rates are invariant under cycHc 
permutations of {A, B, C} , most of the equations below will 
be written for a species a e {A,B, C}, with b and c denoting 
respectively the next and previous species. When q scales as 

q = e-P'^ , 

the dynamics become diffusive: in particular, the probability 
that site 1 < fc < L is occupied by a particle of type a behaves 



Fvo[skit) ^ a] ^ Paik/L,t/L^) , 

where the macroscopic density profiles Paix, r) follow a local 
Fourier's law : if ja{x, r) is the current associated to pa, then 



ja = -dxPa + PPa{Pc ~ Pb) , 

which, together with the conservation law 

drPa = -dxja , 



(4) 



(5) 



gives the hydrodynamic equations [SJl satisfied by the den- 
sity; 



drPa = d^pa + PO^PaiPb - Pc) ■ 



(6) 



These equations conserve the fact that X]a/'a(a^%T) = 1 
(each site is occupied by one of the three species) and 
/p dxpa{x,T) = Va, whcrc ra is the total density of parti- 
cles of type a. The deterministic equations (|4| and ^ are 
only valid in the large L limit, for diffusive time scales, i.e. 

For small (3, the constant density profiles Pa{x) = Va are 
a stable stationnary solution of (|6]l. These constant profiles 
become linearly unstable above a critical value /3„ ll34ll given 
by 



/3* 



-^ with A = 1 - 2 V r^ . 



(7) 



so that the long-time limit of ^ becomes a function aJx) of 
the space variable x for (3 > j3^. It has been argued 1391 (and 
checked numerically) that, in the steady state, these modulated 
profiles do not move. If the phase transition to the modulated 
phase is second order, then it should occur at /? = /?* given 
by (l7]i. A first-order transition may however occur at /3 < /3* 
: this should certainly y^] be the case at least for A < 0, with 



A = E^^-2E^^ 



(8) 



In the following, we study the integrated current QA{t) of 
A particles through the system during time t: 



1 



mW = ^ / dt' [ dx jAix,t')dt' 



"'0 

L I dr I dx JA{x,T)dT . 

Jo 
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FIG. 1 : First cumulant of the integrated current of A particles in the 
ABC model with ta ^ rs = 1/4, for < /J < 2/?.. Measure- 
ments of (Qa) in Monte-Carlo simulations of systems of 80 to 640 
sites are compared to the prediction of the hydrodynamic equation 



This space average fluctuates with time, and we will be in- 
terested in its cumulants, as in Q. Because the difference 
between the space average and the flux through a section re- 
mains bounded, the cumulants of the flux through any section 
are the same as those of this space average in the long time 
limitEl. 



n. MEAN CURRENT 

As the steady-state profiles Pa{x) are time independent ll39ll . 
(|5]l implies that the steady-state currents are constant, ja{x) = 
Ja- They are given by 

Ja = /3 / dxpa{x){pc{x) - Pb{x)) . (10) 



Then, from (|9]l, 

jQAJt)) 13 

t ^ L 



dxpA{x){pc{x)) - psix)). (II) 



{QA{t)) can thus be obtained by calculating numerically the 
long-time limit of ^ and then integrating ( fTOl i. In Figure [T] 
we compare the results of this calculation to numerical mea- 
surements obtained by simulating finite systems of 80 to 640 
sites, for r^i = r^ = 1/4 and < /3 < 2/3*, with /3* given by 
0. 

For /?</?*, the stability of the flat density profiles Pa{x) ~ 
ra leads to a current {QA{t)) ~ |j/3r^(rc —rs)', on the other 
hand, the dependence in f3 becomes non-trivial for /3 > (3^,, 
with a cusp at /3^ . 

For /3 4, /3*, the steady-state profiles are known (see IMI or 
(l27b and (|28l) below) to take the form 

Paix) = ra + v//3-/S*(A'ae2-- + cc.) + 0{f3 - (3,) (12) 
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FIG. 2: Second cumulant of the integrated current of A pailicles in 
the ABC model with ta = tb = 1/4, for < /3 < 2/3,. Mea- 
surements of {Q\)c in systems with 80 to 320 sites are compared to 
those obtained from the macroscopic fluctuation theory l ll3ll4b . 



with known constants Ka, Kb and Kc, leading to an analytic 
expression for {QA{t)) around /3+: 



T<«-(^^f.,(^- 



TB) 



A2 



From this formulation, the generating function of the current 
Q^(i) is the solution of the optimization problem 



log(e 



AQA(t) 



maxL 

Pa ,ja 



t/L 

dr 



dx 



Aj. 



E 



4pa 



(14) 

Finding the density and current profiles Pa{x^ t) and ja{x, r) 
which maximize ( fT4l i is not an easy task. We assume that, for 
large t, the optimum in ( fT4l i is achieved by profiles of fixed 
shape which may drift with a constant velocity v. In order 



^>^QA{t)\ 



to order 2 in 



to obtain (Q^(i))c, we compute log( 
(A, v), before optimizing over v. one then gets optimization 
equations satisfied by these moving profiles, which we solved 
numerically in the case r^ = r^ = 1/4 for < /3 < 2/3*. 

For /3 < /3*, the constant profiles Pa{x) = Ta are still opti- 
mal for A 7^ 0, yielding 



{QAit))c 



2t 



rAil - ta) ■ 



For 13 > (3^,, one has to take into account the dependence of 
the optimal profiles in v and A. This leads (as will be shown 
in the longer 14111 ) to a second cumulant which diverges at 
j3 = p^,. This divergence can be computed exactly thanks to 
the knownledge of the steady-state profiles (fTSl i for /3 | /3*, 
leading to 



(QiW). 
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III. FLUCTUATION THEORY 

The hydrodynamic equations (|6]l describe the determinis- 
tic evolution of the density profiles pa{x,T) in the large L 
limit. For a large, but finite system, one has to take into ac- 
count stochastic corrections. This can be done using fluctuat- 
ing hydrodynamics, where the expression ^ of the current is 
replaced with ifllTgl] 



ja ^qa + -y=7]a{x,T) 



(13) 



where^a = — 9a;Pa+/3pa(Pc— Pfc) istheright-handsideof Q, 
and where the rja are Gaussian noises such that ^^ i]a = 
and 

{ria{x,T)'qa'{x' ,t')) = (Taa'{x,T)6{x -x')5[t -t') , 

with (Jaa = 2/5g(l - Pa) and Gaa' = -'^PaPa' for Q ^ a'. Al- 
ternatively, the macroscopic fluctuation theory ll40ll expresses 
( fT3] l as a large deviation principle, with the probability of ob- 
serving time-dependent density profiles Pa{x, r) given by 



Pro[pa{x,T)] cx cxp 



dxdr \ 



'ipa 



IV. CRITICAL REGIME FOR THE DETERMINISTIC 
HYDRODYNAMICS 

In this section, we analyse how the deterministic equations 
(|6]l, exact in the L ^- oo limit, behave in the neighborhood 
of /3* . The stability analysis of the constant profiles pa = Va 
shows that, as /3 crosses /?* , only the first Fourier modes of the 
Pa become unstable. For /3 close to /?*, one therefore expects 
this mode to relax more slowly than the other Fourier modes. 

One can write, from (|6|l, an evolution equation for these 
slow modes when (3 is close to /?*, and show that they de- 
cay as a power law (instead of an exponential) in the critical 
regime. To do so, we separate in pa{x,T) the first Fourier 
mode, Ra{T)e'^'^^^ , from the other modes, Pa{x, t), so that 



Paix, t) =ra+ Ra{T)e^'''^'' + CC. + Pa{x, t) 



(15) 



ia{x, t) = J air) + —R^e^'-^'' + CC. - dr / dypa{y, t) 
27r Jo 

with Pa '^ Ra '^ i"a- We also suppose that pa{x,T) varies 
slowly, so that drRa = Ra <^ Ra and drPa <C Pa', finally, 
we set 



/3 = /3*(l + 7) 
with 7^1. The leading order of (|4|i then becomes, when 



projected on the first and second Fourier modes, 

2iTTRa =P4Rairc - ^b) +r,(i?, - Rb)] (16) 

dxPa ^MPa{rc - n) + raiPc - Pb) 



+ Ra{Rc-Rb)e 



4i-KX 



CC. 



(17) 



Equation ( fTSI l relates the leading orders of Rb and Rq to Ra, 
so that 

iRB{r)^'-r^^^RA{r)+XBir) 
li?c(r) = 2l^^i±^i?A(r)+xc(r) 

with xb,xc ^ Ra- Equation (fTTI l shows that, at leading 
order, pa is of the form 



Paix,T) = ipae^"'' + CC. 



(19) 



with 



Ai-mpa = P*[Va{rc - n) +ra{ipc - Vb) + Ra{Rc - Rb)] 

whose solution is 



fa 



1 - 2ra ^2 
: n„ 



(20) 



Then, the next-to-leading order of the first Fourier mode of (|4]i 
reads 

-—Ra =2iTrjRa + P^ [{Tc ~n- iVA)Xa + ra{Xc - Xb) 

Ztt 



+ faiR: - Rl) + Rlifc ~ fb)] 



(21) 



(with XA = 0). The Xa can be eliminated by multiplying the 
equation over Ra by {re — rB + i\/A), the one over Re by 
VA, the one over Rb by —rA, and by summing, which yields 



, 2A\Ra\'\ 



(22) 



with A given by ([8]). 

For /? = /?* (7 = 0), the first Fourier mode Ra decays as a 
power law instead of an exponential: 



RAir) 



RAiO) 



1 



A2" 



rA 



\Ra{0W 



with an ampHtude which does not depend on the initial condi- 
tion for large r. 



V. CRITICAL BEHAVIOR OF THE MET 

We now return to the noisy equation ( fT3] l and try to obtain a 
noisy version of (l22T i. By analogy with the deterministic case 
above, we suppose that the first Fourier mode Ra of pa varies 
more slowly, but with a larger amplitude than the other modes, 
so that Ra <C Ra and pa ^ Ra in (fTSl l. When replacing 
dU with ( fTST l ,the leading-order equation ( fT6b is not modified. 



so that (fTsT i still holds : however, the next-to-leading order 
equation (ISTT i is replaced with 

--Ra =^2iTl-lRa + /3* [(re -n- i^/^)Xa + ra{Xc - Xb) 

Ztt 

+ fa{K - Rl) + R*a{fc - fb)] + ^ , (23) 

where the fa are projections of the rja {x, t) on the first Fourier 
mode: 



i^air) = / dxe ^'^''''riaix.T), 
Jo 

so that {Va{T)vl,{T')) ^ CJaa'Sij-T') and {Va{T)Va'{T')) 

0. The second Fourier modes ipa <0 satisfy the equations 



47r 



<Pa = - 4i7rV3a + /3* ['/'a (^c ^ n) + Taiific ~ fb) 



Ra{Rc — Rb)] 



^'^\r) 

Vl ' 



where the z/i ' are projections of the r]a{x, r) as well : hence, 
they fluctuate around their non-noisy expression ( |20] |. In ( l23T l, 
however, these fluctuations (of amplitude 1/vL) are multi- 
plied by Ra, so that they are of smaller amplitude than the 
noisy term Va{T)/VL : therefore, the ipa can be replaced by 
their expression ( |20l ) in ( |23l ). 

Taking a linear combination of ( |23] | to eliminate the Xa as 
in the deterministic case ( I2TI ). we then obtain 



i?A = 47r2 7 



2A |i?^p 



i?A 



/m(t) 



with PA a linear combination of the Va, which verifies 

{pA{T)nAi'r )> = -^ S{t - t ) . 

Finally, the change of variables 



^^M = V FT /(^) 



AL 



with T = Stt t^tt:^ ;= (24) 



A3/2 VI 

leads to a simple rescaled equation: 

f = (7-|/(r)P)/(f)+MT) 

r- A3/2 /3 - /?, 

and with p such that {p{f)p* (f ')) = (5(f — f '). 

Therefore, a system of size L exhibits a critical regime \/3 — 
/3*| ^ 1/VL in which the density profiles pa{x,T) fluctuate 
as sine waves of period 1, with an amplitude scaling as 1/i^/^ 
on a time scale of order XjsfL. The rescaled fluctuations. 
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FIG. 3: Rescaled deviation of the first cumulant of A particles, 
{QA(t)}, from its deterministic hydrodynamics prediction l lllb . as 
a function of 7 ( 125b . for rA~rB = 1/4. Numerical measurements 
for systems of 80 to 320 sites are compared with our analysis of the 
critical point l |27|28t . 



/(t), follow a (complex) damped Langevin dynamics in the 
quartic potential 



V{f) 



-7- 



\.f\' 



4 



and the probability distribution of /(t), P{r, 9, r) = 
Pro[/(r) ~ re*^], satisfies the Fokker-Planck equation 



dfP = -dr 

r 



(r^ -^)r'^P+-rdrP 



^9'oP. (26) 



VI. CRITICAL FLUCTUATIONS OF THE CURRENT 

Let us now discuss the consequences of the slow fluctua- 
tions of the first Fourier mode of the density described above 
on the integrated particle current of a particles, Qa{t)- From 
(fljT l and (flST l, we can express the average instantaneous cur- 
rent, Ja, in terms of /(f) : 



Ja{T) ^ I3ra{n - re) + 2(3{n - r^) 






where Ga, the space average of the noise 7]a, is such that 

{Ga{T)Ga{T')) = 2ra{l - ra)6{T - t') + 0{1/VL) . 



Therefore, the contributions of the fluctuations of the first 
Fourier mode / and of the noise Ga to Ja are of compara- 
ble amplitude. The fluctuations of /, however, occur on the 



slower time scale r ^ t jsTL: hence, they become dominant 
in the integrated current (|9]l 



Qa{t)^L 



t/L^ 



Ja{T)dT , 



with the n-points time correlation function of /(f) giving rise 
to an anormal growth of the n-th cumulant of Qa, {Qa{t))c- 
More precisely, we find that 



(27) 
and 



t STT^y^SAranrc 



with 7 as defined in ( |25] l and with 

1 



A3/2(rfc-rc) 



27rA 



C„(7) 



C„(7) = lim 



dfi..df„(|/(fi)../(f„)|^ 



a time integral of the n-point correlation function of /(f). 

Because /(f) follows a Langevin equation in a quartic po- 
tential (|25] |. one does not have simple analytical expressions 
for the Gnici) for n > 2. Ci(7), on the other hand, only 
depends on the stationary average of |/(f)p: since the sta- 
tionnary state of the Fokker-Planck equation ( |26] ) is P(r, ff) ex 



a-^ir' 



"'' , we obtain easily that 



1 



Ci(7)=7+7T 



2 j7 



dz 



(28) 



Thus, we obtain an analytic expression for the deviation 
of {QA{t)) from its hydrodynamics prediction (fTTT i in the 
neighborhood of /3* . In Figure |3] we compare this expres- 
sion to numerical simulations of systems of 80 to 320 sites, 
for rA = tb = 1/4. 



VII. CONCLUSION 

In this letter we have seen that a driven diffusive system 
like the ABC model, at a second order phase transition, may 
exhibit anomalous Fourier's law at least for the second and 
higher cumulants. This is reminiscent of the cumulants of 
the current which diverge with the system size in the TASEP 
(the totally asymmetric exclusion process) along the first or- 
der transition line 1 121 . The mechanism is however different. 
Here the large fluctuations can be understood by analysing the 
dynamics of the first Fourier mode which becomes unstable at 
the transition whereas, in the TASEP, the large fluctuations of 
the current are due to the presence of the shock. 

The anomalous current fluctuations of the ABC model 
at the phase transition are accompanied by anomalous long 
range density fluctuations which can also be understood in 
terms of the slow noisy dynamics of the first Fourier mode 



(the density fluctuations will be discussed in the forthcoming 
longer version of the present letter l4111 ). 

An interesting open question would be to compare the 
anomalous density and the current fluctuations of the ABC 
model at the transition with those of momentum conserv- 



ing mechanical models, in particular through the dynamics of 
their slow modes. Another interesting question would be to 
study the current fluctuations through other lattice gases (such 
as an Ising model) when there is coexistence of several phases 
at equilibrium. 
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